Experimental analysis of lateral impact on planar brittle material 
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The fragmentation of alumina and glass plates due to lateral impact is studied. A few hundred 
plates have been fragmented at different impact velocities and the produced fragments are ana- 
lyzed. The method employed in this work allows one to investigate some geometrical properties of 
the fragments, besides the traditional size distribution usually studied in former experiments. We 
found that, although both materials exhibit qualitative similar fragment size distribution function, 
their geometrical properties appear to be quite different. A schematic model for two-dimensional 
fragmentation is also presented and its predictions are compared to our experimental results. The 
comparison suggests that the analysis of the fragments' geometrical properties constitutes a more 
stringent test of the theoretical models' assumptions than the size distribution. 

PACS numbers: 46.50.+a, 62.20.M- 



I. INTRODUCTION 

The fragmentation process is a common phenomenon 
that is found both in natural [J, 0, 0, 0, 0, 0, 0] and in- 
dustrial processes || 0, OjJ EH ; and takes place on scales 
ranging from the collision of galaxies [l[ and asteroids 
to the breakup of heavy nuclei @. The underlying 
physics of the fragmentation phenomenon is very differ- 
ent for macroscopic and microscopic systems, since quan- 
tum effects have to be considered in the latter, and also 
because the mechanisms that lead to their breakup are 
not the same. In this work, we focus on the macroscopic 
systems. From both the economic and academic points of 
view, it is important to understand the mechanisms that 
govern this process and, based on them, build models 
that allow predictions for different scenarios. 

As discussed in Ref. brittle materials such as 

glasses and ceramics exhibit no macroscopic plastic de- 
formation as they are subject to tensile loading capable of 
producing a stress level greater than their limiting elas- 
tic limit. Thus, once this stress threshold is reached, 
the fragmentation process immediately starts. The first 
systematic studies of brittle fracture using statistical ar- 
guments were performed about 60 years ago by Weibull 
|13j ]. who introduced a probability distribution function 
(named Weibull distribution) to describe the fragment 
size distributions. 

Recently, this subject has regained a great deal of at- 
tention since Oddershcde et al. [Tj| showed that the size 
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(s) distribution of the fragments produced in the breakup 
of brittle objects is given by a power law s~@ which is 
fairly independent of the specific material employed. It 
is essentially determined by the morphology of the frag- 
menting object and the power law exponent [3 = 1.63 was 
found for spherical gypsum balls, (3 = 1.08 for thin gyp- 
sum disks, and (J = 1.05 for gypsum rods. These authors 
showed that this scaling law is observed for fragments 
whose dimensions are smaller than the smallest dimen- 
sion of the ori gina l object Subsequent investigations 
made in Ref. [la ], where glass rods were dropped from 
increasing heights, suggest that the power law exponent 
rises from (3 rs 1.2 to 1.5 as the violence of the impact 
increases (larger heig hts). A similar experiment has been 
reported in Ref. where a sandwich of thin glass or 
plaster plates, inserted between stainless steel plates, was 
hit by an iron projectile, which was dropped on the tar- 
get at normal angles to its surface. The fragment dis- 
tribution is also given by a power law whose exponent 
increases from j3 = 1.5 to 1.7 according to the violence of 
the impact. Although these studies agree qualitatively, 
the observed power law exponents are somewhat differ- 
ent. Qualitatively different results were obtained in Ref. 
[TtJ where the mass distribution of fragments produced 
in the fragmentation of thick plates of dry clay were de- 
scribed by two power laws of exponents and (3s = 1.5—1.7 
for small fragments and (3l = 1.1 — 1.2 for the larger ones. 
These authors suggest that the exponents arc associated 
with the dimensionality of the fragmenting object and the 
length of the fragments. They also investigated the prob- 
lem using models based on very different pictures for the 
fragmentation process and found that both lead to the 
same mass distributions. This finding strongly suggests 
that the study of other observables, besides the mass/size 
distribution, is necessary to unveil the underlying physics 
of the fragmentation process. This conclusion is one of 
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the motivations for the present work. 

Other recent experimental studies investigated the 
properties of fragmentation in different scenarios. For 
instance, the breakup of closed thin shells due to impact 
and also as a result of the explosion of combustible mix- 
ture has been studied in Ref. (H)]. For practical reasons, 
eggs have been used in these experiments. The authors 
found that the mass distribution follows a power law, in 
both cases, whose exponent (3 = 1.35 ± 0.02 lies between 
Pl and f3s observed in Ref. and mentioned above. 
On the other hand, studies on the fragmentation of mer- 
cury drops [ijj [20| , due to the fall from a fixed height 
on a hard surface, show that the fragment size distribu- 
tion is fairly well described by a power law of exponent 
f3 ss 1 . 1 over a wide range of sizes [l9| . 

The role played by the constraints imposed to the frag- 
menting system has been investigated in Ref. [2l| , where 
slender brittle rods, made up of dry spaghetti (besides 
other brittle materials), have been kept fixed at one end 
and axially impacted at the other. The corresponding 
mass distribution exhibits bumps around A/2 and A/4, 
where A is the preferred wavelength for the buckling in- 
stability. This issue was, to some extent, also investigated 
in Ref. [l6| where the paster plates were laterally bom- 
barded by a high velocity projectile and (3 was found to lie 
within 1.1 and 1.3, which contrasts with the values men- 
tioned above when the plates are hit on their surfaces. 
There have also been experimental studies that focused 
on the internal details of the fractures, such as those re- 
ported in Ref. [22j , where the statistical distribution func- 
tion for the height fluctuations along the fracture length 
was carefully examined and found to be Gaussian. 

Owing to the great complexity of the fracture process 
[23l . I24L [25I . l26l ] , several schematic models have been pro- 
posed to describe the phenomeno n ( see fl9l [HI, US [H, 
M, H, H, HO, HI, [H, Hi, H3, H, HI SlI and references 
therein). Some of them try to incorporate microscopic 
information [25|, US [2?], IH, [37| whereas the majority of 
the approaches are minimalistic models which uses as few 
parameters as possible to describe the fragmentation pro- 
cess [H, [H, [30, HH [H, [H, H3, HE . Due to the very 
small time scale of the fracture dynamics (23|, most of 
the experimental information is usually restricted to the 
final state properties of the system. Therefore, it is very 
difficult to single out the appropriate scenario from these 
models since most of them make predictions which agree 
reasonably well with the available experimental data. 

In this work we present the results of an experiment 
from which we extract further observables associated 
with the fragments' geometry and size, besides their mass 
distribution. By providing more detailed information on 
the fragments produced in the breakup of brittle mate- 
rial, we intend to make it possible to distinguish between 
the different scenarios assumed in the theoretical models 
or, at least, to constrain the range of their free parame- 
ters. We study the fragmentation due to lateral impacts 
on thin plates of alumina and glass, so that our experi- 
ment is similar to one of the measurements made in Ref. 
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FIG. 1: Schematic illustration of the apparatus used in this 
work to study the fragmentation of brittle plates. For details, 
see the text. 

[ill, mentioned above. As described below, we have been 
particularly careful in controlling the impact velocities, 
in order to minimize effects associated with mixing the 
energy deposited into the system. Great care has also 
been taken in order to ensure that the plates are always 
hit along their lateral side, in order to minimize any bias 
due to angle mixing. 

The remainder of this paper is organized as follows. 
Our experimental setup is described in detail in Sect. 
IITI and the model devised to interpret the experimental 
results is presented in Sect. IIIII The results are discussed 
inc Sect. [TV] and the main conclusions are summarized in 
Sect. |Vl 



II. THE EXPERIMENTAL SETUP 

We have used 232 square plates (10cm x 10cm) to 
study the properties of the breakup of thin brittle plates 
due to lateral impacts. Two different materials have 
been considered: alumina and glass. The samples have 
different microscopic structure, but both have noncrys- 
talline atomic arrangement. The alumina targets were 
manufactured by Coorstek Inc. and the glass plates by 
Goes Vidros. None of the plates has been previously syn- 
terized. Their physical properties, as well as the number 
of analyzed plates, are given in Table HI 



TABLE I: Properties of the brittle plates used in this work. 



Material 


Thickness 


Width 


Density 


Quantity 




(mm) 


(mm) 


(g/cm 3 ) 




Alumina 


0.5 


100.0 


0.380 


156 


Glass 


1.0 


100.0 


0.246 


76 



A plate is laid down on a flat surface and it is hit later- 
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FIG. 2: Relation between the applied pressure and the veloc- 
ity of the collision block. 



ally by a piston. No constraints arc imposed on the plate, 
so that it can be freely scattered. A schematic illustra- 
tion of the apparatus is shown in Fig. [1] It consists of 
a compressed air cylinder which is connected to a reduc- 
tion pressure valve. The latter allows one to suitably fix 
the pressure in the pneumatic cylinder connected at its 
end. In order to ensure uniform impacts along the plate's 
length, a 12cmx5cmx2cm (length x height x width) steel 
block is attached to the external end of the piston rod. 
Since the length of this block (12cm) is larger than that 
of the plates (10cm) the impact is very homogeneous. 
The basis for the targets is fixed on an iron plate. To 
it is attached a polystyrene shock absorber bar, which is 
intended to stop the pneumatic cylinder before it reaches 
the piston stroke, thus preserving the physical integrity 
of the apparatus. The compressed air is injected into 
the pneumatic cylinder through the activation of a high- 
speed electrical valve, allowing the steel block to acceler- 
ate and reach the target plate, before being stopped by 
the polystyrene stop. 

The velocity of the block impinging on the target is 
measured using two photodiodes located 1.2cm apart 
from each other. Each of them is illuminated by a laser. 
As the light on either of the photodiodes is interrupted 
by the passage of the block, signals are transmitted to 
a digital oscilloscope, thus allowing the measurement of 
the block speed as a function of the applied pressure. 
The calibration curve obtained with the apparatus is pre- 
sented in Fig. [5J 

The minimal velocity necessary to break the samples 
depends on the material used as target. For our exper- 
imental setup, the threshold for the alumina plates is 
approximately 8.1 m/s and it is 6.4 m/s for the glass 
ones. Below these values, some of the plates occasionally 
break, but produce very few fragments. If a plate does 
not break, it is not used again as it may have suffered 
internal damages. 

Since the targets are not fixed, they move freely af- 



ter having been hit by the piston. In order to confine 
the produced fragments and at the same time preventing 
secondary fragmentation due to further impacts with the 
surroundings, we have placed the system inside a soft 
plastic enclosure with a silk bag at its right end. This 
part of the apparatus is not depicted in Fig. [1] 

The produced fragments are placed on a high resolu- 
tion scanner. This provides detailed images from which 
one determines their geometric properties. We have been 
particularly careful in placing the fragments on the scan- 
ner so as to prevent them from touching each other. For 
the alumina plates, we found that black and white im- 
ages of 600 dpi resolution are suitable for our purposes 
as they allow one to identify fragments whose sizes are 
on the order of ^4 cut = 0.18 mm 2 , i.e. A cut = 1.80 x 10 -5 
smaller than the original objects (Aq). The glass plates 
were scanned as 256 color bitmaps at 200 dpi resolution. 
Therefore, the smallest fragment's area studied in this 
case corresponds to A cut = 1.61 mm 2 . The fragments 
are identified by simply counting the contiguous active 
pixels on the scanned image, similarly to what is done 
in standard cluster recognition algorithms used in per- 
colation theory. We stress that the analysis performed 
in this work provides much more information than the 
traditional ones which focus on the fragments' masses. 



III. A SCHEMATIC MODEL FOR PLAQUE 
FRAGMENTATION 



Before presenting the fragmentation data obtained 
with the apparatus described above, we introduce a 
schematic fragmentation model which will be of help in 
the interpretation of our experimental results. 

As mentioned in the previous section, many models 
have been developed to describe fra g mentation p henom- 
ena in different scenarios [lj, [H, iH, [H, [13, HI, H, HO, 
l3ll . [321 . [33l . [341 . l35l Hil . The confrontation of these models 
with the experimental results presented below is impor- 
tant in order to clarify the essential physical ingredients 
involved in the fragmentation process. However, due to 
the large number of models available in the literature, 
this is beyond the scope of the present work. We therefore 
devised a schematic two-dimensional Monte Carlo model, 
depicted in Fig. [3l inspired in the beautiful experiment 
reported by Xu and Rosakis [231 ] . which we apply to in- 
terpret our experimental results. Although we keep it as 
simple as possible, we introduce the main ingredients we 
consider relevant to the process: 

(i) At the beginning of the process, cracks are created 
at one of the lateral sides of the square plates of 
area Aq — 1 and their propagation directions are 
randomly selected. This number N c of initial cracks 
is closely related to the violence of the impact and 
is a parameter of our model; 

(ii) The propagation of the cracks start simultaneously 
and any of them stops only if its course is inter- 




FIG. 3: Illustration of the propagation of the cracks in the 
framework of the schematic fragmentation model introduced 
in this work. The initial cracks start at the left side of the 
plate. The definition of the bifurcation angle 9 is also depicted 
in this figure. For details, see the text. 
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FIG. 4: (Color online) Fragment multiplicity as a function of 
the initial cylinder pressure. 



glass plates. Then, the probability of creating a new 
crack P c is the only parameter that we vary according 
to the violence of the impact and the material of the 
plates. Owing to the stochastic nature of the model, ten 
thousand events arc run for each case studied below. 



IV. RESULTS 



rupted by another crack or if one of the borders is 
reached. 

(iii) A number Nf of flaw regions are randomly placed 
over the surface of the plate. The number Nf is 
sampled from a Poisson distribution with an average 
value (Nf), which is a parameter of the model. The 
flaw regions are considered to be circles, all of radius 
R, which is another parameter of our model; 

(iv) When a crack enters into a flaw region, a new branch 
is created with probability P c . Its propagation di- 
rection is uniformly sampled between —n/2 and 
+7r/2, with respect to that of the initial crack. The 
latter, continues its course; 

(v) In order to minimize the number of free parameters 
of the model, we argue intuitively that the number 
of initial cracks N c is closely related to the violence 
of the impact, as well as P c . Although it is diffi- 
cult to determine the relationship between N c and 
P c , any reasonable function which gives N c = for 
P c — > and leads to large values of N c for P c — ► 1, 
should provide equivalent qualitative results. For 
simplicity, we adopt N c = — 101n(l — P c ). 

The average number of flaws (Nf) is a parameter of the 
model which is associated with the brittleness of the ma- 
terial. Therefore, different values of (Nf) and R could 
be used for distinct materials. For simplicity, we fixed 
(N f ) = 10000 and R = 0.0005 for both alumina and 



The average fragment multiplicity is displayed in Fig. 
H as a function of the pressure for both alumina and 
glass plates. For each event, fragments of area smaller 
than A cut are not considered as discussed in Sect. [II] For 
consistency, the same procedure is employed in the the- 
oretical calculations. The model results, depicted by the 
upside-down triangles, arc slightly shifted to the right 
in order not to overlap with the data. All the error bars 
shown in this figure correspond to the standard deviation, 
i.e. they are associated with the width of the distribution. 
The model parameter P c was adjusted for each pressure 
used in the experiment. In the case of the alumina plates, 
we used P c = 0.329, 0.393 and 0.425 for P = 10, 12 and 
14 bar, respectively. For the glass plates, we adopted 
P c = 0.53, 0.67, 0.76 and 0.85 for P = 8, 10, 12 and 
14 bar, also respectively. These results show that the 
average fragment multiplicity increases steadily as a func- 
tion of the pressure. One also sees that it is much higher 
for the glass plates than in the case of the alumina tar- 
gets and that, as anticipated in Sect. [TTJ the fragmen- 
tation threshold is much lower in the former case. The 
large variance values, represented by the error bars, re- 
veal that, for a given pressure, the fluctuations are fairly 
large. One should note that, although P c has been ad- 
justed to reproduce the average multiplicity, the model 
also correctly predicts the large variances observed ex- 
perimentally. This is an intrinsic property of the model. 

We now turn to the area distribution and show in Fig. 
El for the alumina plates, the values of F(A), defined as 




FIG. 5: (Color online) Area distribution for alumina plates at FIG. 6: (Color online) Same as Fig. [5] for glass plates. The 
different pressures. The dashed line is a power law of exponent dashed line corresponds to a power law whose exponent is 
P = 1.1. For details, see the text. j3 = 1.2. For details, see the text. 



F(A) = jJ™n(A>)dA>, (1) 

where n(A)dA is the number of fragments with area be- 
tween A and A + dA. 

The dashed line shown in this figure represents the 
power law fit F(A) oc A~P . The results show that the 
experimental data can be fairly well approximated by 
this function over about three decades with (3 = 1.1. 
This exponent is in agr eement with some of the results 
reported previously [la , [l^ | . One observes a steeper drop 
of F(A) at large areas, which becomes more important 
as the violence of the impact increases. Nevertheless, 
the fragment distribution associated with not too large 
fragments (A/A < 10~ 2 ) remains essentially unchanged 
within the pressure range studied in this work. 

The agreement of the model results, depicted in this 
figure by the full lines, with the experimental data is very 
good. Small discrepancies are observed only for large ar- 
eas, but they are compatible with the experimental un- 
certainties in this area region. 

Qualitatively similar results are observed for the glass 
plates, as is shown in Fig. [S] The exponent of the power 
law function used in this case is slightly larger than that 
obtained above and it corresponds to (3 = 1.2. The range 
over which this function remains a good approximation 
to the actual behavior observed experimentally is con- 
siderably smaller than in the case of the alumina plates. 



The suppression of large areas is much more pronounced 
in the present case. This is in qualitative agreement with 
the results shown in Fig. [¥J since the fragment multiplic- 
ity is much larger in the case of the glass plates than for 
the alumina ones. 

The model also reproduces the area distribution very 
well in this case. However, as mentioned in Sect. U it 
was demonstrated in Ref. [l7[ that different assumptions 
for the fragmentation mechanisms may lead to equiva- 
lent mass distributions. Therefore, the description of 
F(A) should be regarded as a preliminary selection crite- 
rion and further comparisons with the experimental data 
should be made in order to validate any model. In the 
following, we examine some size correlations in order to 
seek for further information on the fragmentation pro- 
cess. 

The average value of the largest area, (^.Largest), within 
each event is displayed in Fig. [7] as a function of the frag- 
ment multiplicity for the alumina and glass plates. In 
both cases, (^Largest) decreases as a function of the frag- 
ment multiplicity for not too high multiplicities. In the 
case of the alumina plates, one observes a slight increase 
at the highest multiplicity values. This clearly contrasts 
with the behavior observed in the glass plates. We be- 
lieve that this deviation is due to the poor statistics of 
the very high multiplicity events since no anomalies were 
observed in these events. Nevertheless, it could also be 
explained by a geometrical preference for the creation of 
the fragments. Since this information is not available in 
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FIG. 7: (Color online) Average value of the largest area within 
each event as a function of the multiplicity. In the case of the 
alumina plates, multiplicity bins of 11 units were used whereas 
21 units were employed for the glass plates. 




Model 



the present measurements, this hypothesis cannot be in- 
vestigated at this moment. The predictions of the model 
agree fairly well with the alumina data for multiplici- 
ties smaller than 70 - 80. It decreases smoothly as a 
function of the fragment multiplicity. In the case of the 
glass plates, small discrepancies are observed at the low- 
est multiplicities, where the model tends to overpredict 
the size of the largest fragment whereas it underprcdicts 
this observable for the alumina plates. This is an indi- 
cation that, in spite of the very good overall agreement 
observed above, the details associated with the fragment 
distribution might not be correctly reproduced by the 
model. 

Important information on the fragmentation pattern 
may be obtained through the Dalitz plot [111, which is 
based on the properties of the largest fragments. More 
specifically, this is constructed by calculating 

3 

Xi = Ai/Y,Ak, k= 1,2,3, (2) 

k=l 

where {A^} corresponds to the area of the largest three 
fragments within each event. The quantity \i repre- 
sents the perpendicular distance to the k-th side of an 
equilateral triangle of height 1, into which a point asso- 
ciated with a given event is plotted. By construction, 
all the points lie inside the triangle. Due to geometri- 
cal constraints, an event point is univocally defined by 
the pair of distances (xi,Xj)i * 3- The indices {k} 
are randomized in each event in order to eliminate ar- 
tificial structures, i.e., the labels assigned to each of the 
three largest fragments are shuffled and, therefore, we ex- 
plicitly avoid systematic correlations associated with the 
fragment sizes. A large concentration of points close to 
the vertices indicates a fragmentation mode in which the 



FIG. 8: (Color online) Dalitz plot of the area distribution 
associated with the fragments produced in the fragmentation 
of the alumina plates. For details, see the text. 



size distribution has a big fragment whereas the others 
are appreciably smaller. Three big fragments of approx- 
imately the same size lead to points grouped around the 
center of the triangle. On the other hand, when two frag- 
ments have approximately the same size, but are much 
larger than the third one, one finds points gathered close 
to the middle point of the triangle's sides. 

The experimental Dalitz plot displayed in Fig. [5J ob- 
tained in the fragmentation of the alumina plates, show 
that there is an important contribution from events in 
which one fragment is much larger than the others. The 
results shown in this figure correspond to the cylinder 
pressure equal to 14 bar, but the conclusion remains valid 
for the lower pressures employed in this work. The cor- 
responding model results are also displayed in the lower 
part of this figure and the qualitative agreement with the 
experimental features is, once more, very good. 

A different behavior is observed in the case of the glass 
plates, whose Dalitz plot tends to give points which are 
grouped near the center of the triangle, for all the pres- 
sures we considered. As discussed above, this means that 
the three largest fragments have approximately the same 
area. Figure [5] shows the corresponding experimental and 
theoretical plots. Due to the low experimental statistics 
for this kind of plot, the data associated with different 
pressures are grouped in this figure, which does not affect 
our analysis as they are very similar for all pressures. The 
model correctly predicts the tendency to suppress contri- 
butions associated with a dominant fragment. Therefore, 
the agreement with the data is also fairly good in this 
case. 

Further insight into the properties of the fragment dis- 
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FIG. 9: (Color online) Same as fig.|S]for the glass plates. For 
details, see the text. 



tribution can be obtained by examining the shape pa- 
rameter Q, that we define as: 
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FIG. 10: (Color online) Average Q value versus the fragments' 
area for different values of the pressure for the fragmentation 
of the alumina plates. For details, see the text. 



where P denotes the perimeter of the fragment and A 
corresponds to its area. Large Q values indicate elon- 
gated shapes whereas approximately circular or squared 
shapes give Q close to unity. In this way, this quan- 
tity allows one to obtain information on the shape of the 
fragments. It is worth mentioning that we analyzed the 
fractal dimension of the borders and found that it is es- 
sentially one dimensional. Figure [TTJ displays the average 
value of Q for different impact velocities, obtained in the 
fragmentation of the alumina plates. It shows that the 
small fragments tend to be fairly elongated but this ten- 
dency is quickly weakened as the violence of the impact 
increases. On the other hand, the not too small frag- 
ments, i.e. A/Ao > 10~ 4 , are much less elongated than 
the small ones, for all the pressures studied in this work. 
The understanding of this property requires the devel- 
opment of a model that describes the fracture process 
accurately. In this context, it should also be interesting 
to investigate whether there exists a preferential direc- 
tion to the elongation, but this is beyond the scope of 
this study as the present experimental setup does not 
provide this kind of information. 

The average value of Q obtained with our model is 
also shown in Fig. [10] and is depicted by the upside-down 
triangles. The agreement with the experimental values is 
fairly good for A/Aq > 10 -3 , but the model predicts too 
elongated fragments for smaller areas. We come back to 



this point below. 

The fragments produced in the breakup of the glass 
plates have different properties since the (Q) versus A 
curve is fairly fiat, except at the lowest pressures, as 
shown in Fig. [TTJ It should be noticed that the verti- 
cal scales in Figs. [TU] and [TTJ are not the same. More 
specifically, the elongation of the fragments originated 
from the glass plates is much smaller than that observed 
in the case of the alumina objects. Despite the exper- 
imental uncertainties, it is clear that the model, once 
more, systematically predicts too elongated fragments as 
already noted in the case of the alumina plates. These re- 
sults suggest that the observable Q might be a useful tool 
to selecting models which give an appropriate picture to 
the fragmentation process. 

The tendency to form too elongated fragments in our 
model can be appreciably reduced by not allowing the 
cracks initiated at the flaw points to propagate in di- 
rections close to the original crack. We have checked 
that the quality of the agreement with all the observ- 
ables presented above is maintained if the bifurcation 
angle 9 is, for instance, restricted to 7r/6 <| 9 |< tt/2. 
All the other parameters keep their values, except for 
P c which is slightly changed (by less than 4%) in some 
cases. To illustrate this fact quantitatively, the Q values 
obtained in this case are depicted by the squares in Figs. 
[TOl and [TTJ It is clear that the agreement with the data 
improved considerably. Since there are experimental and 
theoretical studies [2J, [2f| [39|, H(| that give support to 
this angular restriction, our results strongly suggest that 
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FIG. 11: (Color online) Same as Fig. 1101 for the glass plates. 



the new cracks start at relatively large bifurcation angles. 
In spite of these encouraging results, we preferred not to 
tune this model parameter and we keep this conclusion 
on a qualitative level. We postpone detailed discussions 
on this subject to future work when further observables 
will be analyzed. In this work we intend to stress that 
models can reproduce global quantities, such as the mass 
distribution or the average fragment multiplicity, whereas 
they can fail in describing more detailed information on 
the fragments' properties and that important physical 
aspects can be revealed by restricting the model param- 
eters. 

As a final remark, it would be fair to speculate whether 
the experimental (Q) values, for small areas, are not bi- 
ased by the fact that we discard fragments whose areas 
are smaller than A cut . It should be noticed that this 
lower limit is a safe cut which is much larger than the 
image resolution. The specific values of A cut have been 
selected for the alumina and glass plates based on the 



criterion that the dimensions of the fragments arc larger 
than the thickness of the plates due to the ambiguity 
associated with the identification of the corresponding 
dimensions. Although smaller fragments could be identi- 
fied, they have been excluded from the actual calculation 
of (Q) and other observables. Qualitative analyses of 
such small fragments indicate that our results should not 
be impacted by the consideration of these fragments. 
V. CONCLUDING REMARKS 

In summary, we have presented the results of an exper- 
iment in which the fragmentation of two different brittle 
materials due to lateral impact was studied. In agree- 
ment with former experimental works, we found that the 
size distribution is given by a power law over a wide range 
of sizes. The power law exponents are very similar for 
both materials, i.e. glass and alumina, and respectively 
correspond to j3 R* 1.2 and and y 1.1, which are close 
to the values obtained in Refs. [3, EH, E3, EH ■ 

Going beyond those studies, the present experiment 
also focused on the geometric properties of the fragments. 
The description of these properties provide strict tests to 
theoretical models. More specifically, we found that the 
size distribution and average fragment multiplicities are 
very well described by the simple fragmentation model 
presented here, whereas it failed to reproduce the aver- 
age elongation of the fragments' shapes. This quantity 
is reproduced by the model only if the bifurcation an- 
gle of the fractures is restricted to relatively large angles 
(tt/6 <| 8 |< 7r/2) with respect to the propagation of the 
initial crack. This is in qualitative a gree ment with the 
experimental results reported in Ref . 12311 a nd the theo- 
retical studies discussed in Refs. j2ol l39l |40| . We believe 
that careful comparisons of these experimental observ- 
ables [4l| with the predictions of different models might 
be very useful in helping to establish a clear scenario for 
the fragmentation process. 
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